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ABSTRACT 

We study gravitational instability and consequent star formation in a wide range of isolated disk 
galaxies, using three-dimensional, smoothed particle hydrodynamics simulations at resolution suffi- 
cient to fully resolve gravitational collapse. Stellar feedback is represented by an isothermal equation 
of state. Absorbing sink particles are inserted in dynamically bound, converging regions with number 
density n > 10 3 cm~ 3 to directly measure the mass of gravitationally collapsing gas available for star 
formation. Our models quantitatively reproduce not only the observed Schmidt law, but also the 
observed star formation threshold in disk galaxies. Our results suggest that the dominant physical 
mechanism determining the star formation rate is just the strength of gravitational instability, with 
feedback primarily functioning to maintain a roughly constant effective sound speed. 
Subject headings: galaxy: evolution — galaxy: spiral — galaxy: star clusters — stars: formation 



1. INTRODUCTION 

Stars fo rm in galaxies at hugely varying rates 
Il998ajh The mechanisms that control the 

;as are widely de- 
2002: 



( Kcnm 

star for mation rate from interstellar 
bated ll Shu. Adams k Lizanoi Il98fl: lElmeereenl ._ ._. 
iLarsonl 120031 iMac Low k Klessenl |2004|) . Gravitational 
collapse is opposed by gas pressure, supersonic turbu- 
lence, magnetic fields, and rotational shear. Gas pres- 
sure in turn is regulated by radiative cooling and stel- 
lar and turbulent heating. Despite this complexity, 
star-forming spiral galaxies follow two empirical laws. 
First , stars only form above a critical gas surface den- 
sity ijMartin k KennicuttH2001[) that appears to be de- 
termined bv the lTo"mnreTi|l964j) criterion for gravitational 
instability. Second, the rate of star formation i s propor- 
tional to a power of the total gas surface density ijSchmidtJ 
Il959t IKennicuttlll998bl) . 

A number of groups have simulated disk galaxies in 
isolat ion or in mergers , or in cosmological c o ntexts , 
e.g.. iMihos k HernauistJ <1994D: iFriedli Bend" Jl995[k 
Sommer-Larsen. Gelato k Vedell lll999MSpringe]l <l200Cfl: 
Barnesl (|2002lk iGovernato et all l)2004() . iRobertson et all 
( 2004) review this work. However, in these simulations, 



star formation is generally set up with empirical recipes 
a priori. The origin of the observed Schmidt law remains 
unclear. 

Recent cos mological simula tions with moderate mass 
resolution bv lKravtsovl (|200Hfl show that the Schmidt law 
is a manifestation of the overall density distribution of 
the ISM, and find little contribution from feedback. How- 
ever, the strength of gravitational instability was not di- 
rectly measured in his work, so a direct connection could 
not be made between instability and the Schmidt Law, as 
we do here. The importance of gravitational instability 
in controlling larg e-scale star formation was emphasized 
bv lFriedli k Bend l)1995j) and lElmegreenl lj2002j) . A simi- 
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lar conclusion comes from the observation that thin dust 
lanes in galaxies only form in gravitationally u nstable 
regions t| Dalcanton. Yoachim k Bernstei n 2004). 

We simulate a large set of isolated galaxies to inves- 
tigate gravitational instability and consequent star for- 
mation. In this Letter, we examine star formation as 
a function of gravitational instability, and compare the 
global Schmidt law and star formation thresholds derived 
from our simulations to the observations. 

2. COMPUTATIONAL METHOD 

We use the smoothe d particle hydrodynamics 
fSPH ) code GADGET ij Snrmgel. Yoshida k White! 
120011). modified to include absorbing sink particles 
ijBate. Bonncll k Pried 1 19951) to directly measure the 
mass of gravitationally collapsing gas. Sink particles, 
representing star clusters (SCs), replace gravitationally 
bound regions of converging flow that reach number 
densities n > 10 3 cm~ 3 . (These regions have pressures 
P/k ~ 10 7 K cm -3 typical of star-forming regions.) 

Our galaxy model consists of a dark matter halo, and 
a disk of stars and isothermal gas. The galaxy structure 
is bas ed on the analytical work by I Mo. Mao k W hite 

( 1998) , as i mplemen t ed nu merically bv lSpringel fc W hite 

(1999) and lSpringeil l)2000|) . The isothermal sound speed 
is chosen to be either c\ = 6 km s _1 in models with 
low temperature T or C2 = 15 km s _1 in high T models. 
Table 1 lists the most important model parameters. The 
Toomre criterion for gravitational instabil ity that couple s 
stars and gas, Q sg is calculated following Rafikov ( 2001), 
and the minimum value is derived using the wavenumber 
fc of greatest instability and lowest Q sg at each radius. 

The gas, halo and disk particles are distributed with 
number ratio N g : Nh : Yd = 5 : 3 : 2. The gravita- 
tional softening lengths of the halo eh = 0.4 kpc and 
disk ed = 0.1 kpc, while that of the gas e g is given 
in Table 1 for each model. The minimum spatial and 
rklessej^iR-^olutions in the gas are given by e g and twice 
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TABLE 1 

Galaxy Models and Numerical Parameters 



Model a 


k h 


Q S g(LTT 


Qs g (HT) d 






m s s 


G50-1 


1 


1.22 


1.45 


1.0 


10 


0.08 


G50-2 


2.5 


0.94 


1.53 


1.0 


10 


0.21 


G50-3 


1.5 


0.65 


1.52 


1.0 


10 


0.37 


G50-4 


9 


0.33 


0.82 


1.0 


10 


0.75 


G100-1 


1 


1.08 




6.4 


7 


0.10 


G100-1 


1 




1.27 


1.0 


10 


0.66 


G100-2 


2.5 




1.07 


1.0 


10 


1.65 


G100-3 


4.5 




0.82 


1.0 


10 


2.97 


G100-4 


9 




0.42 


1.0 


20 


5.94 


G120-3 


1.5 




0.68 


1.0 


20 


5.17 


G120-4 


9 




0.35 


1.0 


30 


10.3 


G160-1 


1 




1.34 


1.0 


20 


2.72 


G160-2 


2.5 




0.89 


1.0 


20 


6.80 


G160-3 


4.5 




0.52 


1.0 


30 


12.2 


G160-4 


9 




0.26 


1.5 


40 


16.3 


G220-1 


1 


0.65 




6.4 


15 


1.11 


G220-1 


1 




1.11 


1.0 


20 


7.07 


G220-2 


2.5 




0.66 


1.2 


30 


14.8 


G220-3 


4.5 




0.38 


2.0 


40 


15.9 


G220-4 


9 




0.19 


1.0 


10 


16.0 



a First number is rotational velocity in km s 1 at virial radius. 
b Percentage of total halo mass in gas. 

c Minimum initial Q sg for low T model. Missing data indicates 
models not run at full resolution. 
d Minimum initial Q sg for high T model 
°Millions of particles in high resolution runs. 
'Gravitational softening length of gas in pc. 
g Gas particle mass in units of 1O 4 M0. 



the kernel mass (~ 80m g ). We adopt typical values for 
the halo concentration parameter c = 5, spin parame- 
ter A = 0.05, and Hu bble constant Hq = 70 km s _1 
Mpc^ 1 l|Springe]| 12000). The spin parameter used is a 
typical one for galaxies su bject to the t idal fo rces of the 
cosmological background. jReed et all (2003) suggest a 
wide range of c for galaxy-size halos. However, this pa- 
rameter is based on a simple model of th e halo forma- 
tion time ( Navarro. Frenk fc White! 119971). with noorlv 
known distribution <|Mo et al. Hl998j)Tspringel fc White! 
(1999) suggest that c = 5 is theoretically expected for 
fiat, low-density universes. 

Models of gravitational collapse must satisfy three 
numerical criter i a: t he Jeans resolutio n criterion 
(iBate fc BurkertJ 119971 hereafter BB97; iWhitwortbl 
1998), the gravity-hydro balance criterion for gravi- 
tational softening (BB97), and the equipartition cr i- 
terion for par t icle m asses l|Steinmetz fc White! Il997|) . 
iTruelove et alJ l)1998l) suggest that a Jeans mass must 
be resolved with far more than the Nk — 2 smoothing 
kernels proposed by BB97. Therefore we performed a 
resolution study of model G100-1 (LT) with N tot = 10 5 , 
8 x 10 5 , and 6.4 x 10 6 , corresponding to iVk ~ 0.4, 3.0 and 
23.9, respectively. We find convergence to within 10% of 
the global amount of mass accreted by sink particles be- 
tween the two highest resolutions, suggesting that the 
BB97 criterion is sufficient for the problem considered 
here. 

We performed 24 simulations satisfying all three crite- 
ria, including six models of low mass galaxies with low T 
to study the effect of changing the effective sound speed. 
We also set a minimum value of N to t > 10 6 particles for 



lower mass galaxies resolved with fewer particles. 

3. GLOBAL SCHMIDT LAW 

To derive the Schmidt law, we aver age Ssfr and £ 
over the star forming region following iKennicutn |l989), 
with radius chosen to encircle 80% of the mass in sinks. 
To estimate the star formation rate, we make the as- 
sumption that individual sinks represent dense molecular 
clo uds that form sta r s at s ome efficiency. Observations 
by iRownd fc Youna ((1999) suggest that the local star 
formation efficienc y (SFE) in molecul ar clouds remains 
roughly constant. iKennicutil l)1998bf) shows a median 
SFE of 30% in starburst galaxies dominated by molecu- 
lar gas. This suggests the local SFE of dense molecular 
clouds around 30%. We therefore adopt a fixed local SFE 
of e = 30% to convert the mass of sinks to stars. Note 
that this local efficiency is different from the global star 
formation efficiency in galaxies, which measures the frac- 
tion of the total gas turn ed into stars. The global SFE 
can range from 1-100% l|Kennicuttlll998bj) . depending 
on the gas distribution and the molecular gas fraction. 




(Mo pc" 2 ) 



Fig. 1. — Schmidt law from fully resolved low (open symbols) 
and high (filled symbols) T models listed in Table HI that showed 
gravitational collapse. The colors indicate the galaxy rotational 
velocities, while the symbol shapes indicate the gas fractions, as 
specified in t he legend. The bla ck line is the best fit to the obser- 
vations from Kcnnicutt 1 1998b), while the red line is the best fit 
to the simulations. 



Figure ^ shows the Schmidt law derived fro m our sim- 
ulation s. The best fit to the observations bv iKennicuttJ 
(1998b) gives a Schmidt law Ssfr = ^4Sg as with global 
efficiency A = (2.5±0.7) x 10 -4 and power law a — 1.4± 
0.15, where Ssfr ^ s g iven m units of Mq kpc~ 2 yr _1 , 
and S gas is given in units of M Q pc -2 . A least-squares 
fit to the models listed in Table ^ (both low T and high 
T) gives A = (1.4 ± 0.4) x 10~ 4 and a = 1.45 ± 0.07, 
agreeing with the observations to within the errors. 

Note LT models tend to have slightly higher SF rates 
than equivalent HT models. Thus, observations may 
be able to directly measure the effective sound speed 
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(roughly equivalent to velocity dispersion) of the star- 
forming gas in galactic disks and nuclei. More simula- 
tions will be needed to demonstrate this quantitatively. 

Our chosen models do not populate the lowest and 
highest star formation rates observed. Interacting galax- 
ies can produce very unstable disks and trigger vigorous 
starbursts (e.g., I'LL Mac Low fc Klessedl2004|) . Quies- 
cent normal galaxies form stars at a rate below our mass 
resolution limit. Our most stable models indeed show no 
star formation in the first few billion years. 

4. STAR FORMATION THRESHOLD 

A threshold is clearly visible in the spatial distribution 
of gas and stars in our galaxy models, as illustrated in 
Figure El The critical value of the instability parameter 
at threshold can be quantitatively measured from the ra- 
dial profile as indicated in the middle panel, which shows 
a sharp drop of Esfr at i? ~ 2i?d- The critical values of 
Q sg and Q g at the threshold i? t h are shown in the bot- 
tom panel of Figure |3 for all the fully resolved models 
listed in Tabled The critical values of <5 sg appear to be 
generally higher than Q g in the same galaxy, and both 
have lower values (< 1) in more unstable models. 

Most galaxies not classified as starbursts have gas frac- 
tions comparable to or less than our most stable mod- 
els, so the observation of a threshold value of Q g ~ 1.4 
may reflect the stability o f the galaxies in the sample 
l)Martin fc KennicuttH200lt) . Observed variations in the 
threshold also appear to occur naturally. If we only use 
the Toomre criterion for the gas Q g we get slightly larger 
scatter than if we include the stars and use the combined 
criterion Q sg , but the effect is small. 

5. DISCUSSION AND SUMMARY 

What controls star formation in different galaxies? 
Our models suggest the answer is the nonlinear devel- 
opment of gravitational instability. Figure [3] shows the 
correlation between the star formation timescale tsf and 
the initial minimum Q sg (min) for fully resolved models 
listed in Table □ The best fit is t sf = (34 ± 7 Myr) x 
exp [(4.2 ± 0.3)Q sg (min)]. Quiescent star formation oc- 
curs where Q sg is large, while vigorous starbursts occur 
where Q sg is small. Thi s differs from the emphasis on su- 
personic turbulence bv IKravtsovl lj2003T) . The maximum 
strength of instability Q sg (min) depends on the mass of 
the galaxy and the gas fraction. The larger the halo 
mass, or the larger the gas fraction, the smaller resulting 
Q S g (mm), and thus the shorter rsp. 

Typical observed starbu rst times of 10 8 y r are consis- 
tent with our fit for t sf ilKennicuttlll998bl). This also 
agrees with the observations bv [Mac Art hur et alJ (2004) 
that the star formation rate depends on the galaxy poten- 
tial. iMcGaueh et all ((2000) show a break in the Tully- 
Fisher relation for galaxies with V c < 90 km s _1 , sug- 
gesting a transition at this scale. Indeed, our models 
with V c < 100 km s _1 and gas fraction < 50% of the 
disk mass appear to be rather stable (Q sg > 1.0), with 
no star formation in the first 3 Gyrs, while models with 
V c > 120 km s^ 1 become less stable, forming stars easily. 
This is also consistent with the rotatio nal velocity above 
which dust lanes are observed to form ijDalcanton et al.1 
2004). 

We have deferred inclusion of explicit feedback, mag- 
netic fields, and gas recycling to future work. However, 
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Fig. 2. — (Top) Star formation threshold illustrated by the low 
T model G220-1 with TVtot = 6.4 X 10 6 . Log of gas surface density 
is shown, with values given by the color bar. Yellow dots indicate 
SCs, while the red circle shows Rth- (Middle) Radial profiles of 
star formation rate (yellow circles), and Toomre Q parameters for 
stars Q s (asterisks), gas Q g (circles), and stars and gas combined 
Q sg (diamonds). The red line shows i? t h- Bottom: critical values 
of Q sg (filled symbols) and Q g (open symbols) at H t h for both low 
(red) and high (black) T models. 



we believe each will have minor effects on the ques- 
tions considered here. The assumption of an isother- 
mal equation of state for the gas implies substantial 
feedback to maintain the effective temperature of the 
gas against radiative cooling and turbulent dissipation. 
Real interstellar gas has a wide range of temperatures. 
However, the rms velocity disper sion generally falls 
withi n the range 6-12 km s _1 (e.g.. lElmegreen fc Scald 
2004). Direct feedback from starbursts may play only 
a minor role in quenching subsequent star formation 
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Fig. 3. — Star formation timescale rgp as a function of ini- 
tial Q sg (rain), for both low T (open symbols) and high T (filled 
symbols) models. The solid line is the least-square fit. 



(e.g.. IKravtsovl I2003t IMonacd |200l) . perhaps because 
most energy is deposited not in t he disk but above 
it as superbubbles blow out ( e.g.. iFuiita et all |200^t 
lAvillez fc Breitschwerdll2004|L iKim fc Ostrikerl i|2001fl 

demonstrate that swing and magneto- Jeans instabilities 
operating in a gaseous disk occur at Q g ~ 1.4, suggesting 
that magnetostatic support is unimportant. The lack of 
gas recycling both from disrupted molecular clouds and 
from massive stars will change the detailed patterns of 
star formation, but probably not the overall results. 

Simulations o f iso lated, isothermal disks by 
IRobertson et al.l l|2004() show large-scale collapse in 
their centers leading to disks far smaller than observed, 
which they argued was caused by an isothermal equation 
of state. This behavior does not occur in our model 
with physical parameters close to theirs, but resolu- 
tion sufficient to resolve the Jeans length. Similarly, 



iGovernato argue that several long-standing 

problems in galaxy simulations such as the angular mo- 
mentum catastrophe may well be caused by inadequate 
resolution, or violation of the other numerical criteria. 
We will present more resolution studies in future work. 

In summary, our models reproduce quantitatively not 
only the Schmidt law, but also the star formation thresh- 
old in disk galaxies. We find a direct correlation between 
the star formation rate and the strength of gravitational 
instability. This suggests that gravitational instability in 
effectively isothermal gas may be the dominant physical 
mechanism that controls the rate and location of star 
formation in galaxies. Unstable galaxies were more com- 
mon at early cosmic times, so our results, together with 
merger-induced starbursts (iLi et al.| |20041 m ay account 
for the Butcher-Oemler effect l|Butcher fc Oemlerlll984f) 
of increasing blueness of galaxies with redshift. Massive 
galaxies form stars quickly, which may account for the 
downsizing effect that star formation first occurs in big 
galaxies at hig h redshift, while modern starburst galax- 
ies are small JCowie et alJ Il996t iPoggianti et alJl2004t 
iFerreras et al.ll2004|) . The slow evolution of star forma- 
tion in our low mass models rese mbles that observed 
in low surface brightness galaxies l)van den Hoek et alJ 
2000). 
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